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1. Introduction 

In recent years, active control of sound and vibration in aeroelastic panel 
has gained a great deal of attention due to many possible applications to 
aerospace and related industries. In the absence of a flow field, such problems 
have been studied by several authors. For an overview of this subject, one 
is referred to the paper [1] — [4]. Most recent advances in active control of 
sound and vibration can be found in the conference proceeding [5]. 

_ In the current research project, we seek to investigate some problems in 
active control of panel vibration excited by a boundary- layer flow over a flat 
plate. In the first phase of this investigation, we have studied the optimal 
control problem of vibrating elastic panel induced by a fluid dynamical load- 
ing. For a simply-supported rectangular plate, the vibration control problem 
can be analyzed by a modal analysis. The control objective is to minimize 
the total cost functional, which is the sum of a vibrational energy and the 
control cost. By means of the modal (eigenfunction) expansion, the dynam- 
ical equation for the plate and the cost functional are reduced to a system 
of ordinary differential equations and the cost functions for the modes. For 
the linear elastic plate, the modes become uncoupled. The control of each 
modal amplitude reduces to the so-called linear regulator problem in control 
theory. Such problem can then be solved by the method of adjoint state. 
This method was used successfully in our previous study on the control of 
thermal fluctuation in a shear flow [6]. The optimality system of equations 
was solved numerically by a shooting method. The results of this investi- 
gation will be summarized in this report. We have also begun to look into 
the control problem for a nonlinear panel vibration. Analytically we have 


1 



obtained some preliminary results. Numerical solution of the optimal con- 
trol of a nonlinear panel (or beam) in one space dimension is in progress. 
Further development of the vibration control problem will be detailed in the 
next report. 

2. Results on Active Control of Panel Vibrations 

Let us consider a viscous flow past over the elastic panel. The flow is 
governed by the well know Navier-Stokes equation: 

^ + («. V)u = --Vp + i/Au, (1) 

at p 

where the notations are standard. For a slightly compressible flow, the con- 
tinuity equation reads 

j- t + V • (pu) = 0. (2) 

The panel is regarded as an elastic plate with thickness /i, Young’s modulus E 
and Poisson’s ratio 7. Under a uniform tension with N > 0 (or compression 
with N < 0) and the fluctuating wall pressure, the vertical displacement ( 
of the plate satisfies the following equation: 

^^ = JVAC-DA 2 < + p„+ ? (i,() (3) 

where p w is the plate density, p w the wall pressure fluctuation, and q is the 
applied force as the active control. The constant D is the stiffness of the 
plate defined by 


12(1 — 7 2 ) 
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According to the boundary-layer theory, given an upstream velocity field 

U, the flow near the plate can be determined by the Prandtl’s approximation. 

In particular, if the panel is located on the x — y plane, the pressure gradient 
dt) 

w ~~~ across the boundary-layer is nearly constant, where x = ( x,y , z). Suppose 
Oz 

that the mean flow outside the boundary-layer is parallel to the plate so that 
U = (f/oo, 0,0) + Ui(x,t), where U\ is a small perturbation. To derive the 
equations for acoustic quantities u,pi and p\ , we let 

u = u 0 + Ux , p = Po + Pi and p = p 0 + p x (5) 

where u 0 ,Po and p 0 are flow variables associated with the mean flow. As in the 
stability analysis, we introduce a parallel flow approximation. Then, in view 
of (5), one obtains the acoustic equations from (1) and (2) by linearization: 

^7 + (uo • V)U! = — “Vpi, (6) 

dt po 

+ V • {po^i + piUo) = 0. (7) 

For an isentropic flow, pi and p\ are related by 

Pi = Pi/c\ (8) 

where c is the speed of sound for the unperturbed flow. Aside from a static 
displacement, the vibration of the panel is described by the perturbation w 
of equation (3) as follows: 

r \ 2 

= NAw - DA 2 w + f + q(x,t), (9) 
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where f = p w is the fluctuating part of the wall pressure excited by the 
unsteady boundary-layer flow. The coupling of the acoustic equations (6) and 
( 7 ), and the plate equation is through the boundary conditions. For the plate 
equation (9), since it is simply supported by a periodic structure, we need 
only to analyze the problem over a fundamental domain 0 < x < a, 0 < y < b 
and impose the boundary conditons: 

w(x,y,t) = 0 at x = 0, a; y = 0,6. (10) 

dp 

Since the pressure gradient = 0 across the boundary, the wall pressure 
can be determined from the perturbed potential flow field U i through an 
approximate Euler’s equation, that is, 

P„ = F(U l ). ( 11 ) 

To counter this excitation, a control force q(x,y,t) was introduced in ( 9 ). 
The objective of the active control is to minimize the average vibrational 
energy and the control cost: 

J(q) = -L F [ M^) 2 + fl(Aw ) 2 + 7 |VH 2 + kq 2 }dtdxdy, (12) 

where the time T may be infinite, D is the basic domain {0 < x < a ,0 < 
V < b}; a,/?, 7 and k are positive constants. In the language of the optimal 
control of a distributed parameter system, the equation ( 9 ) is known as the 
equation of state and J(q) defined by (12), the objective or cost functional. 
Here the physical problem of vibrational control reduces to an optimiza- 
tion problem: Given the wall pressure excitation p w , find an optimal control 
q*{x, y, t ) from a certain admissible class Q of functions which minimizes the 
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objective functional J(q ), that is, 

J(q*) = min{ J{q), q in Q}. (13) 

To obtain an analytical solution, we consider the case of simply supported 
boundary conditions: 

w(x,y,t ) = 0 at x = 0,a; y = 0,6, 


= 0 at x = 0,a and i,t) = 0 at y — 0, b. (14) 

ox* oy 


The initial conditions are given by 


dw 


w(x, y, 0) = w 0 (x, y), -^-(x, °) = w i( x ’ »)■ 

It is well known that the set of functions 


(15) 


. . , m7r . mr 

<p mn (x, y) = 2 sin x sin —y, m,n = 1, 2, • • • 

a b 


(16) 


are orthogonal eigenfunctions associated with the plate equation (9) and the 
corresponding eigenvalues are 


A mn = N[(^y + (^) 2 ] + D[(^y + (^)T- 


,nic s 


,m7i\ 

T 


(17) 


In terms of the above eigenfunctions, we can expand the displacement w, the 
wall pressure / and the control as follows: 


m,n= 1 


(18) 


/(#,$/,£)— 3/) 

m,n= 1 


(19) 
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and 


oo 

q{ x ,yj)= J2 9m, n{t)<Pm n {x,y), 

m,n= 1 

where the coefficients w mn etc. are computed by 


w r 


m — ( ) — f [ ^(*£? t/, (^? t/)d£<f?/ , 

JO Jo 


(20) 


and so on. A substitution of the expansions (18)-(20) into the equations (9), 
(15) and (12) yields the following uncoupled system of equations: 


and 


{ Pw^mn “h ^mn^mn — fmniX) + 9ran(0 ? 

^mn(O) = ^0,mn 5 ^mn ( 9 ) = ^1 ,mn 


( 21 ) 


OO 

J(q) = E -M?) , (22) 

m,n=l 

where 


Jnn(9) = l + fyiLtO}** 4 ( 23 ) 

for ra,n = 1,2, • * Since the modes are uncoupled, if the cost J mn for each 
mode is minimized, so does the total cost J. 

For a given (m, n)-mode, dropping all the subscripts, we are led to con- 
sider the so-called “linear regulator” problem in optimal control: Find the 
control q in the equation 

pw + Xw = f(t ) + q(t) , 

1 . 

iu(0) = w 0 , iu(0) = u>i , 
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which minimizes 


J(q) = jjj; jf {“•“ 2 + (' u|2 + kq 2 )dt , 


(25) 


where // mn is given as in A mn with D and N replaced by /? and 7, respectively. 
By the method of adjoint state,^, for the cost to be minimal, the state w 
and its adjoint v must satisfy the optimality system: 


t 

pw + \w = f(t) + l(aw - V ), 
w(0) = u>(0), th(0) = wi, 


(26) 


and 


{ pv + Xv = (aA + p)™, 
v(T) = v(T) = 0. 


(27) 


The optimal control q* is given by q — — n). One notes that, due 

to the coupling between w and v, the above system (26)-(27) is a two- 
point boundary- value problem. Numerically it can be solved by the shooting 
method. Some numerical results for the original modal equations (21) and 
(23) have been obtained. 

For example, we choose a = b = w, tu 0 = W\ = 0 and T = 4, and set 


= 7 T \ — 77 cos(m 2 -I- n 2 ) 1/2 t, m,n = 1,2, 

(m 2 + n?) 

All the parameters are taken to be one except for /?, which is zero. The 
maximal amplitude of vibration under a optimal control has been computed 
and some results, corresponding to 4 modes (m + n = 4), are shown in Fig. 
1 to Fig. 3. In the above figures, the solid curves represent the controlled 
amplitudes, which are in contrast with the uncontrolled ones. It is seen that 
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the control is very effective in reducing the vibration amplitudes. For an 
independent interest, the controlled mode shape at t = 4 is plotted as shown 
in Fig. 4. 


3. Active Control of Panel Fluttering 

The problem of panel flutter is a classical problem in aeroelasticity. The 
problem is of great concern to those who are involved in the design and 
operation of airplanes. The fluttering is a phenomenon of nonlinear self- 
excited panel vibrations. To include the nonlinear effect, we shall confine 
our analysis to the case for which the panel is long in the ^-direction so 
that the problem involves only one space-dimension x. Here in contrast with 
Eq. (9), the equation of motion for the panel reads: 
d 2 w 


dt 2 


= (N + N x )^-D^ + f(x,t) + q(x,t), 


dx 2 


dx 4 


(28) 


w = 


w(x, 0 ) 


d 2 w 

-Q -2 =0 at X = 0,a, 

dw ( n\ 

W °’ -df( X ’°' = Wl ’ 


where 


2 a Jo \dx ) 

is the nonlinear effect due to streching, E and h are the Young’s modulus 
and the panel thickness, respectively. Without the active control q , the 
nonlinear panel vibration due to the aerodynamic forcing / has been studied 
extensively. However, a comprehensive study of the active control problem 
is still lacking. Subject to the optimality criterion J(q) given by Eq. (12), we 
have studied this problem analytically by a modal expansion. To this end, 
we first derived the optimality system for the state equation (28), in which 
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the optimal control must satisfy the following terminal-value problem: 


where 


d 2 q 

} dP 


, .. d 2 q ,d 2 u> ,d 2 w 

(N + N x )-?— — D— - — 2 8{ , q)- 


dx 2 dx 4 ' dx 2 

1 . d 2 w r,d*w d 2 w 

+ TA a ~fo2 +7 all) > 


dx 2 


'ax 2 


9 = 


8 2 q 

= 0 at x = 0,a, 

0 and = 0 at t = T, 

ax fc/9 at 


8 = (— — ) and 
v 2a ' 

(w,q) = I wqdx. 

Jo 


( 29 ) 


Since the optimality system (28) and (29) is now nonlinear, the modal ex- 
pansion (18) - (20) is x alone yields an infinite system of coupled equations. 
We approximate the system by a modal truncation. The truncated system is 
sovled numerically. The numerical study is still in progress. The numerical 
solution will be presented in the next report. 
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